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ABSTRACT 

Hot Jupiters are giant planets on orbits a few hundredths of an AU. They do not share their system 
with low-mass close-in planets, despite these latter being exceedingly common. Two migration chan¬ 
nels for hot Jupiters have been proposed: through a protoplanetary gas disc or by tidal circularisation 
of highly-eccentric planets. We show that highly-eccentric giant planets that will become hot Jupiters 
clear out any low-mass inner planets in the system, explaining the observed lack of such companions to 
hot Jupiters. A less common outcome of the interaction is that the giant planet is ejected by the inner 
planets. Furthermore, the interaction can implant giant planets on moderately-high eccentricities at 
semimajor axes < 1 AU, a region otherwise hard to populate. Our work supports the hypothesis that 
most hot Jupiters reached their current orbits following a phase of high eccentricity, possibly excited 
by other planetary or stellar companions. 

Subject headings: planets and satellites: dynamical evolution and stability — stars: individual (Kepler- 
18, Kepler-23, Kepler-58, Kepler-339) 


1. INTRODUCTION 

Hot Jupi ters were among the f irst exoplanets to be 
discovered (|Mavor fo Quelozlll995ll . However, their ori¬ 
gin is still not understood, and models for their mi¬ 
gration history fall into two categories: “Type II” 
migration at early times through th e protoplanetary 
gas disc dLin et al.1119961 lWardl1l997ll ; and migration 
at late times as planets’ eccentricities are excited by 
gravitational scattering in packed multi-pl anet systems 
(jRasio fc Fordlfl996t iChatteriee et al.l[2008H . or by secu¬ 
lar perturbat ions from more distant planetary or binary 
companions dWu fc Murravll2003t IWu fe Lithwickl [20111 : 

IBeauge fe Nesvornvll2012t lPetrovichll20l4 I2015D . High- 
eccentricity migration may better explain the observed 
misalignments betw e en stellar sp i n and planetary orbits 
dTriaud et, all 1201(11: iWinn et al l I201(H IWu fc Lithwickl 

120111 IBeauge fc Nesvornvl l2012t iStorch et all 120141 as 
well as the innermost semi-majo r axes of the bulk 
of the hot Jupiter population dFord fe Rasiol 120061 
iPlavchan fc Bilinskil 120131 iValsecchi fc Rasiol 1 201 4) 1 A 
requirement of this channel is that hot Jupiters have 
(or had in the past) planetary or stellar companions 
on wide orbits, and indeed recent studies estimate that 
around 70% of hot Jupit ers have companion giant planets 
or st ars on wide orbits (jKnutson et al.112014 iNgo et afl 

EioTl . 

On the other hand, hot Jupiters are not found to 
have low-mass, close-in companions. No such compan¬ 
ions have yet been found by radial-velocity surveys, while 
survey results from the Kepler spacecraft found no ev¬ 
idence of additional transiting companio ns or transit 
timing variations in hot Jupiter systems dSteffen et al.1 
I20T2D : this deficit was statistically significant compared 
to multiplicities of warm Jupiter and hot Neptune sys¬ 
tems. Nor have ground-based searches for companions 
that may ca use strong transit timing variations prove d 
fruitful (e.g.. lHover et alJl2012tlMacieiewski et al.ll20T3h . 


despite these being sensitive to Earth-mass companions 
in mean motion resonance with a hot Jupiter. How¬ 
ever, low-mass planets on close orbits are extremely com¬ 
mon around stars that do not host hot Jupiters: re¬ 
sults from Kepler transit photometry show that 52% 
of stars have at least one planet with P < 85 days 
and I?pi > 0.81?® dFressin et al.1 l2013h : while radial- 
velocity surveys similarly show that 23% of stars host 
at least one pl anet with P < 50 days and m p i > 3M® 
([Howard et al.1 .2010;. Furthermore, such planets of¬ 
ten occur in multiple systems: the statistics of Ke¬ 
pler candidate multiplicities require s a si gnificant contri- 


bution from multi-planet systems 

Lissauer et al.l 201 lk 

Fang & Margotl 120121; iFressin et al.l 

2013ft. In many svs- 


tems, then, migrating giant planets that will become hot 
Jupiters will interact with formed or forming systems of 
low-mass planets. 


The lack of close companions to hot Jupiters can 
help to dist i nguis h the different migration modes 
dSteffen et all I2012T ) . Simulations show that a giant 
planet migrating through an inner gas disc to become 
a hot Jupiter does not ne cessa ril y suppress planet for¬ 
mation in the inner disc (iMandell fc Sigurdsson 1200.4 
iFogg fo Nelsonll2005ll2007a«P. 120091: IMa/ndell et a!.ll2007lh 

while embryos migrating after the giant form a resonant 
chain behind it and may accrete into a planet of de¬ 
tecta ble size (iKetchum et al.1 1201 It lOgihara et al.112013L 
l20lil . 

In contrast, we show in this paper that during high- 
eccentricity migration, the giant planet almost always 
destroys all low-mass planets on orbits of a few tenths 
of an AU. Previous studies have shown that scattering 
among multiple giant planets can clear out material in 
the terrestrial planet region around 1 AU, t hrough di¬ 
rect scattering dVeras fc Armitage 2005, 2006) or secular 
resonance sweeping dMatsumura et, al.ll2013li . and that it 
can suppress terrestrial pla net formation in this region 
(1 Raymond et al. 1 120111 1201 21b We choose to focus our at¬ 
tention on very close-in systems more relevant for com- 
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parison to Kepler observation (~ 0.1 AU), which may 
have significant mass in inner planets (up to ~ 40 Mg in 
total). We further consider the general case of a highly- 
eccentric giant planet, which may represent the outcome 
of scattering but which may also arise through other ec¬ 
centricity excitation mechanisms such as Kozai pertur¬ 
bations or other secular effects. 

In Section 2 of this paper we briefly review the popu¬ 
lation of planet candidates revealed by the Kepler space¬ 
craft. In Section 3 we describe the numerical approach 
we take to study the interaction of eccentric giant plan¬ 
ets with close-in inner planets. In Section 4 we present 
the results of our numerical integrations. We discuss our 
findings in Section 5, and conclude in Section 6. 

2. PLANETARY MULTIPLICITIES 

We show the multiplicities of the population of 
Kepler planet candidates by taking the catalogue 
of Kepler Objects of Interest (KOIs) from the Ql- 
Q16 data release at the NASA Exoplanet Archive 
(NEA) http://exoplanetarchive.ipac.caltech.edu/ (re¬ 
lease of 2014-12-18; accessed 2015-01-08). This provided 
a list of 7348 planet candidates. KOIs may be genuine 
planets or false positives, with false positive probabil- 
ities up to 1 in 3 in some regions of para meter space 
(|Santerne et al.l l2012t iCoughlin et al.1 I2014D . Moreover, 
parameters for some planets in the NEA are unphysical. 
We therefore performed several cuts on this list to at¬ 
tempt to remove false positives and poorly-characterised 
candidates: 

No FPs: Removal of any candidate classed as a false 
positive in the NASA Exoplanet Archive (in either of the 
columns “disposition using Kepler data” or “Exoplanet 
Archive disposition” V 5139 candidates. 

L+ll: Following ILissauer et all ( 201 1). we consider 
only planets with SNR > 16, P < 240 days and 
R < 22.41?0, thus ensuring completeness and removing 
candidates with unphysically large radii. 3678 candi¬ 
dates. 

L+l l & no FPs: Applies the cuts from ILissauer et all 
1 201 l b and also removes any false positives identified in 
the NEA Q1-Q16 data. 3228 candidates. 

NEA good: Removes NEA-identified false positives, 
and furthermore only includes planets that are listed as 
“confirmed” or “candidate” in at least one of the disposi¬ 
tion columns, ensuring that the planets, if not confirmed, 
have passed some vetting to ensure a low probability of 
a false positive. 2052 candidates. 

Plots of planet radius versus period for these four sam¬ 
ples are shown in Fig [TJ where the contrast between soli¬ 
tary hot Jupiters and sociable low-mass planets is ap¬ 
parent. Although the numbers of single versus multiple 
systems vary (in particular, the NEA good sample has 
a great many multiples, as it is heavily influenced by 
the va lidat ion of numerous multiple-candidate systems 
bv lRowe et al.ll2014D . for our purposes the key observa¬ 
tion is that hot Jupiters at the top left of the plot are 
single. We do see some candidate hot Jupiters with com¬ 
panions, but these detections are not robust. For exam¬ 
ple, with the L+ll and L+ll & no FPs cuts, we find 
KOI-199.01 and KOI-199.02. The latter component is 
marked as a background eclipsing binary in the Q1-Q6 
data from the NEA. In the NEA g ood s ample, we find 
KOI-338, confirmed bv lRowe et al.1 (I2014T) as Kepler-141. 


TABLE 1 

Stellar parameters used in our 

INTEGRATIONS 


Name 

M*/Mq 

R*/Rq 

Reference 

Kepler-18 

0.972 

1.108 

i 

Kepler-23 

1.11 

1.52 

2 

Kepler-58 

0.95 

1.03 

3 

Kepler-339 

0.902 

0.802 

4 


Referen c es. — (1) [Cochran ct al. (20TJ); (2) 
(3) IStefFen et all tl 201 31") : (4) 


This object has an unphysically large stellar radius in the 
NEA, measuring 191?0, l arge r than each of its planets’ 
orbital radii. iRowe et all (l2014h assign a radius of 0.8 J?q, 
reducing the radii of the planet candidates proportion¬ 
ately. Hence, we do not consider either of these potential 
exceptions to our assertion that hot Jupiters are single 
to be reliable. 

We adopt the L+ll & no FPs sample as the most 
reliable. This has 3228 planet candidates, forming 
2136 single-planet systems, 282 doubles, 109 triples, 31 
quadruples, 13 quintuples and 2 sextuples. We restrict 
attention to the triples and lower multiplicites as they 
offer better statistics. 

3. NUMERICAL METHOD 
3.1. N-body Model 

We conduct an extensive ensemble of N- body integra¬ 
tions with the Mercury package dChambersI 119991 1. We 
consider a highly-eccentric giant planet interacting with 
systems of three low-mass planets at ~ 0.1 AU, chosen 
from among Kepler triple-candidate systems, assuming 
that the three transiting planets are the only ones present 
in the inner system. Our systems are representative of 
the range of planet sizes of the multi-planet Kepler sys¬ 
tems (Fig[TJ). Integrations are run for 1 Myr. 

We adopt the Bulirsch-Stoer algorithm with an error 
tolerance of 10~ 12 . Within the 1 Myr integration dura¬ 
tion, energy conservation is generally good; we reject a 
small number of runs with A E/E > 10 -3 . Collisions 
between bodies are treated as perfect mergers, and we 
consider a planet ejected from the system if it reaches 
a distance of 10 000 AU from the star. The code does 
not incorporate general-relativistic corrections, but this 
is unimportant as the dynamics is dominated by scatter¬ 
ing. 

For our main integration runs, we take three-planet 
systems from the Kepler triples and add to the system 
a highly-eccentric giant planet with a small pericentre. 
Our exemplar Kepler systems are Kepler-18, Kepler-23, 
Kepler-58 and Kepler-339. Kepler-18, -23 and -339 all 
have planets with orbits from ~ 0.05 to ~ 0.12 AU, and 
span the range of planetary radii of the Kepler multi¬ 
ple systems. Kepler-58 has planets on somewhat wider 
orbits, 0.09-0.23 AU. These systems are marked in the 
space of Kepler candidates in Fig [TJ 

The Kepler photometry allows a direct determination 
only of planet radii, but masses are more significant dy¬ 
namically. Where available, we have taken masses de¬ 
termined by transit timing variations or radial veloci¬ 
ties. Where these are unavailable, we have estimated 
masses based on a mass-radius or density-radius rela- 
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Fig. 1.— Periods and radii of candidate Kepler single-, double- and triple-transit systems, with various selection criteria. 

A: No FPs —Known false positiv es removed; B: L+ll —Cuts on radius, SNR and period applied following I Lissauer et alJ <120111 ): C: L+ll 
& no FPs — [Lissauer et al.fd201U) cuts applied and false positives removed; D: NEA good —Only retains those candidates confirmed or with 
favourable vetting. See text for further details. Hot Jupiters are alone in systems, in contrast to both smaller planets and giant planets at 
longer orbital periods. We adopt the L+ll & no FPs selection as the most reliable. We mark the triple-candidate systems selected for our 
N-body integrations in blue, with their Kepler numbers shown. 


tion (|Weiss fe M arcv 2014). System parameters used for 
the simulations are given in Tables |T] and [2] 

For each of these systems, we conducted integration 
suites with different properties of the giant planet. For 
all four systems, we conducted integrations with the gi¬ 
ant planet’s initial semi-major axis of 10 AU, while for 
Kepler-18 and -339 we also conducted integrations start¬ 
ing at 1.25 AU. Within each combination of system and 
semi-major axis, we conducted 21 sets of 256 integra¬ 
tions, one set for each pericentre value q from 0.01 AU 
to 0.20 in steps of 0.01 AU, and a final set at 0.25 AU 
(see Fig ED. Within each set, half of the giants were on 
prograde and half on retrograde orbits; within each sub¬ 
sample, the orientation of the orbit was isotropic in the 
respective hemisphere. The giant was always released 
from apocentre. The giant’s mass and radius were set 
to Jupiter’s values. Our set-up assumes that during the 


initial excitation of the giant planet’s eccentricity, there 
is no effect on the inner system, a reasonable assumption 
(for example, a tightly-packed system o f plan ets pr otects 
itself against the Kozai effect. Unnanen et al.lll997 ). 

For the inner systems, we assigned the planets ini¬ 
tially circular orbits with inclinations of up to 5° from 
the reference plane, giving a maximum mutual inclina- 
tion of 10° with the d istribution peaking at around 3.5° 
(j Johansen et al.ll2012fl . We conducted an additional in¬ 
tegration suite for the Kepler-18 system starting from a 
highly-coplanar configuration of inner planets (inclina¬ 
tions up to 0.001°), finding little impact on the outcome. 
The initial orbital phases of the inner planets were ran¬ 
domised. We conducted some integrations without giant 
planets to verify that our three-planet systems do not 
destabilise themselves on relevant timescales. No unsta¬ 
ble systems were found over lMyr (128 runs for each 
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TABLE 2 

Planetary parameters used in our integrations 


Name 

a/AU 

Mpl/M® 

Rpi / -^® 

Reference 

Kepler-18 b 

0.0477 

6.9 a 

2.00 

i 

Kepler-18 c 

0.0752 

17.3 

5.49 

i 

Kepler-18 d 

0.1172 

16.4 

6.98 

i 

Kepler-23 b 

0.0749 

T86 3 

1.89 

2 

Kepler-23 c 

0.0987 

8.05 a 

3.25 

2 

Kepler-23 d 

0.125 

5.60 a 

2.20 

2 

Kepler-58 b 

0.0909 

18.0 

2.78 

3 

Kepler-58 c 

0.1204 

17.5 

2.86 

3 

Kepler-58 d 

0.2262 

7.33 a 

2.94 

4 

Kepler-339 b 

0.0551 

3776 s 

1.42 

4 

Kepler-339 c 

0.0691 

1.74 a 

1.15 

4 

Kepler-339 d 

0.0910 

1.86 a 

1.17 

4 


Referen ces. — (1) ICochrjmejualJ |2011j) ; (2 j lForr] et al.l 
11201211 : LSI IWu fc L ithwick (2013); (4) Rowe c t al.lJ2^^fi 


a Mass is not measured directly: estimated from a mass—radius 
or density—radius relation JWeiss Sz Marc"vll2014D . 

Kepler triple studied). 

We also conduct some ancillary integrations to test the 
effects of the relative orbital energies of the inner planets 
and the giant planet on the probability of ejecting the 
giant. All of these integrations were performed with a 
semimajor axis of 10 AU and a pericentre of 0.02 AU for 
the giant. We tested two hot Jupiter syste ms (51 Pegasi, 
iMavor fc dueled 119951 iButler et alJ 120001: and HAT-P- 
7. lPal et al.l 120081 1 and one high-mul tiplicity system dis - 
covered by radial velocity (r Ceti; iTuomi et al.ll2013T) . 
We also conducted additional integrations for Kepler-18 
with the giant’s mass set to 0.1, 0.3 and 3 Mj at 10 AU, 
and with the giant’s mass set to 1 Mj at 5 and 2.5 AU. 
The ejection fractions from these integrations are used in 
the discussion of the effects of orbital energy on ejection 
probability (see §4), but their statistics are not otherwise 
discussed. 

3.2. Tidal Model 

Although we do not incorporate tidal forces into our 
N-body integrations, we post-process the planets sur¬ 
viving at the end of the 1 Myr N-body integration to 
follow their orbital evolution under tidal forces. To 
model the tidal evolution of the planets after the inter¬ 
action between the giant and the inner planets has con¬ 
cluded, we use th e sim p le “const a nt Q” model in the 
form given in iDobbs-Dixon et ahl (|2004D . We include 
only the planetary tide, which is the most important for 
the planets’ eccentri city decay until the host star leaves 
the Main Sequence (lYillaver et al.ll2014lh We adopt val¬ 
ues for the planets’ tidal quality factors of Qpj = 10 6 
for the giants; Q' pl = 10 5 for the “Neptunes” Kepler- 
18c, d and their merger products; and Q' pX = 10 2 for 
the super-Earth Kepler-18b. These values are at the 
high end of those estim ated for Solar System giants 
dGoldreich &; Soteil Il96(t) but _comparable to estimates 
for exoplanets (1 Jackson et al.l 1200811 . We tidally evolve 
our systems for 10 Gyr. We note that observed systems 
may have had less time to tidally evolve, and a shorter 
evolution with a proportionally smaller Q will give the 
same outcome. 


4. RESULTS 

Our N-body simulations show that, in most cases, the 
systems resolve to one of two outcomes on time-scales 
much shorter than the integration duration: either all of 
the inner planets are destroyed (usually by collision with 
the star), leaving a single eccentric giant; or the giant is 
ejected by the inner planets, leaving 1 -3 inner planets in 
the system, all of low mass. Examples of orbital evolution 
leading to these outcomes are shown in Fig El 

For our chosen systems, we explore varying the giant 
planet’s pericentre and semimajor axis (Fig [3]). So long 
as the giant’s orbit is intersecting at least one of the 
inner planets’, the majority of integrations lead to one 
of the two outcomes in less than 1 Myr (Fig QJ. For 
tidal circularisation of the giant’s orbit to form a true 
hot Jupiter, the pericentre must be a few hundredths 
of an AU (see below); within this distance, nearly all 
of our simulations result in one of the two outcomes of 
ejection of the giant or destruction of all the inner plan¬ 
ets. The overwhelming outcome is that the three inner 
planets are destroyed, most commonly by collision with 
the star, although in some cases the giant accretes one 
or more, which may significantly enrich the core of the 
giant. Indeed, many hot Jupiters are observed to have 
enriched cores (e.g.. iGuillot et aLll2006D . After evolving 
our surviving giant planets under tidal forces, we find 
that giants that have accreted one or more inner plan¬ 
ets are more likely to become hot Jupiters than those 
that have not, although the relative infrequency of colli¬ 
sion means that most of the hot Jupiters we form have 
not accreted other planets (Table [3]). Due to the extreme 
collision velocities at these small orbital radii, the giant’s 
radius may be inflated by colliding with a smaller planet 
llKetchum et alJ 12011 ). while the impact velocity when 
two smaller planets collide can be several times their es¬ 
cape velocity, meani ng that collisions may generate copi¬ 
ous debris dLeinhardt fc Stewartll2012lh The inclination 
of the giant planet with respect to the inner planets does 
not have a strong effect on the outcome, although with 
a retrograde giant the fraction of destroyed inner planets 
colliding with the giant rather than the star rises slightly, 
as does the number of coexisting systems when the initial 
q < 0.10 AU. 

While the ejection of a Jovian planet by Neptune-sized 
ones may seem surprising, the ratio of ejections of the in¬ 
coming giant to destruction of the inner planets can be 
understood in terms of the orbital energies of the two 
components (Fig[5j): as the orbital energy of the giant 
is decreased (whether through lower mass or through 
higher senrinrajor axis), ejection becomes more likely. 
Planets scattering from near-circular orbits at semimajor 
axes of ~ 0.1 AU would be in a regime favouring colli¬ 
sions, as their ph ysical radius is a significant fraction of 
their Hill radius dJohansen et alJ l2012t iPetrovich et alJ 
120141) : equivalently, the ratio of their escape velocity to 
orbital velocity is small, meaning that orbits are not per¬ 
turbed as much during close encounters. However, for 
the highly-eccentric planets we consider here, ejection is 
easily achieved because a small transfer of energy from 
the inner planets to the giant can lead to a significant 
change in the latter’s semimajor axis. Ejection of the gi¬ 
ant is a common outcome for the most massive systems 
of inner planets we consider when the giant comes in on 
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0 1 xio 4 2X10 4 3x10 4 4x10 4 5x10 4 IxlO 4 2x10 4 3x10 4 

t [yr] t [yr] 

Fig. 2.— Example orbital evolution from our integrations. Semimajor axis (a), pericentre ( q ) and apocentre (r max ) are shown 
for the giant and the three inner planets in two runs from the Kepler-18 integrations. The giant starts with a pericentre of 0.02 AU, 
meaning its orbit is overlapping those of all the inner planets. A: the eccentricities of the inner planets are stirred up, while the giant 
planet’s semimajor axis fluctuates as a result of energy exchange with the inner planets. Just after 40 kyr, the giant planet and planet c 
collide; shortly after that, the remaining inner planets are driven into the star. The outcome is the giant planet surviving alone on a more 
tightly-bound orbit, with little overall change in its pericentre. B: the inner planets’ orbits begin intersecting sooner, and c and d collide 
just after 3 kyr. Around 6 kyr, the giant is kicked onto a very wide orbit, only completing a few more orbits before acquiring a hyperbolic 
orbit at around 27 kyr. Meanwhile the c—d merger product swallows planet b, leaving a single inner planet in the system after ejection of 
the eccentric giant. 


TABLE 3 

Numbers of planets forming hot Jupiters in each 
SYSTEM AFTER 10 GYR OF TIDAL EVOLUTION. 



Form HJ 

Do not form HJ 

Kepler-18, a = 10 AU 

Accreted planet 

Did not accrete planet 

146 (22%) 
150 (13%) 

509 (78%) 

987 (87%) 

Kepler-18, a = 1.25 AU 

Accreted planet 

Did not accrete planet 

204 (36%) 
407 (11%) 

365 (64%) 

3174 (89%) 

Kepler-23, a = 10 AU 

Accreted planet 

Did not accrete planet 

90 (17%) 
322 (12%) 

445 (83%) 

2344 (88%) 

Kepler-58, a = 10 AU 

Accreted planet 

Did not accrete planet 

49 (12%) 
178 (8%) 

361 (88%) 

2108 (92%) 

Kepler-339, a = 10 AU 

Accreted planet 

Did not accrete planet 

138 (26%) 
434 (22%) 

384 (74%) 

1559 (78%) 

Kepler-339, a = 1.25 AU 

Accreted planet 

Did not accrete planet 

184 (40%) 
507 (14%) 

280 (60%) 

3114 (86%) 


Note. — “Hot Jupiter” is here defined to be a planet that 
becomes tidally circularised to e < 0.1. 


a wide orbit, but is rare when the inner planets are less 
massive or the giant’s semimajor axis is smaller. 

Giant planets that destroy the inner planets experi¬ 
ence some change to their orbital elements (Fig [6]) . Peri- 
centres may change slightly, while semimajor axes may 
be significantly reduced. Many of the surviving giants 
maintain the small pericentres needed for tidal circulari- 
sation, and will become hot Jupiters after long-term tidal 
evolution: after 10 Gyr, between 8% and 23% of giants 
in our integrations circularise to e < 0.1, depending on 


TABLE 4 

Weighted numbers of hot and warm Jupiters from our 

SIMULATIONS, AFTER IOGyR OF TIDAL EVOLUTION. 


System 

a < 0.1 AU 

0.1 AU < a < 0.76 AU 

Kepler-18 

10675 

4732 

Kepler-18, a = 1.25 AU 

23108 

8474 

Kepler-23 

14740 

6184 

Kepler-58 

8292 

4024 

Kepler-339 

20273 

6020 

Kepler-339, a = 1.25 AU 

24557 

7605 

KOIs 

189 

125 

RV-detected 

677 

255 

RV-detected, e > 0.4 


80 


Note. — We weight each planet by l/(a(l — e 2 )) to correct for the 
geometric transit probability. Below the line we show the numbers of 
KOIs with R > 8 i?® from the Q1-Q16 KOI list with the L+ll & no 
FPs cuts, as well as a sample of RV-detected giant planets (see text). 

the inner planet configuration and the initial semimajor 
axis of the giant planet (see Table [3]). The final semima¬ 
jor axes of these hot Jupiters are < 0.06 AU, implying 
pre-circularisation pericentre distances (after interaction 
with the inner planets) of < 0.03 AU. 

We also form a population of giant planets with large 
pericentres (q > 0.05 AU, too large for tidal circular- 
isation), relatively small semimajor axis (a < 1 AU), 
and moderately high eccentricity (e > 0.5). It is hard to 
populate this region through in-situ scattering of close-in 
giant planets t hat may have migrat ed through a proto¬ 
planetary disc (|Petrovich et al.lDOMl) . since planets scat¬ 
tering from these semi-major axes are inefficient at ex¬ 
citing eccentricity from circular orbits as they are in 
a regime favouring collisions; see the discussion above 
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Fig. 3.— Outcomes of our simulations, after 1 Myr of orbital evolution, shown as stacked bars as a function of the 
giant planet’s initial pericentre q. Blue shows systems where the giant destroys the inner planets (possibly by accreting one or more 
of them), red/orange shows systems where the giant is lost (usually by ejection, although collisions with the star can occur at q = 0.01 
AU) and greys show systems where the giant and at least one inner planet coexist after 1 Myr. In a small number of systems, all planets 
are lost (yellow), or energy conservation is poor (black). Red triangles mark the initial semimajor axes of the inner planets. Yellow and 
red circles show the systems to scale (planet radii inflated by a factor of 50). When the giant’s pericentre comes inside the semi-major axis 
of the outermost inner planet, most integrations result in ejection of the giant or destruction of the inner planets. Destruction of the inner 
planets is favoured unless the inner planets are very massive (as in Kepler-18 and -58) and the giant is weakly bound (a = 10au). 
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q [au] 


Fig. 4. Time until loss of the first planet in our Kepler- 
18 runs with an initial giant planet semimajor axis of 

10 AU. We overplot the time spent at each pericentre when un¬ 
dergoing Kozai cycles from a 1 Mq perturber at 1000 AU: in many 
cases, planets are lost within the time it would take the giant’s 
pericentre to pass through the relevant region. 

and iPetrovich et al.i (l2014f) . Nor is it straightforward 
to populate this region through tidal circularisation, 
as the planets lie below the tidal circularisation tracks 
down which planets move on astrophysically interesting 
timescales. Planets may enter this region as a result of 
“stalled” Kozai migration, if the parameters of the body 
driving the planet’s eccentricity are sufficient to continue 
driving Kozai cycles during the tidal dissip ation process 
(iDong et al.ll20l4lDawson fc Chiangi2014lh or as a result 
of secular chaos (|Wu fc Lithwickll2011fl . Either of these 
pathways entails certain constraints on the perturber ex¬ 
citing the eccentricity, and in particular it is not clear 
to what extent the conditi ons needed to tr igger secular 
chaos are met in practice (|Davies et al.ll20l1) . Our model 
of a high-eccentricity giant interacting with inner planets 
permits us to populate this same region, without relying 
on suitable parameters of the exciting body. Unfortu¬ 
nately most Kepler candidates do not have measured ec¬ 
centricities, but still we can compare the numbers of giant 
planets in semi-major axis bins: from our simulations, 
after 10 Gyr of tidal evolution we find around 2-3 times 
more giant planets in the range a € (0,0.1) AU than 
in a G (0.1,0.76) AU (corresponding to a 240d period) 
after correcting for the geometrical transit probability 
while in our Kepler sample, we find only around 50% 
more (Table IU . Hence, our results are consistent with 
the observed population, as we might expect the “warm 
Jupiter” region beyond 0. 1 AU to b e populated to some 
extent by disc migration dLin et al.lll996t> — which better 
explains the low-eccentricity warm Jupiters—while some 
of the hot Jupiters may be destroyed a s a res ult of tides 
raised on the star (IValsecchi fc Rasiol 12014) . We can 
also consider the planet population detected by radial- 

1 Ratios are similar if we stop the tidal evolution after 1 Gyr, 
although with Kepler-58 and Kepler-18, a = 10 AU the ratio is a 
little lower at 1.7. 



•giant/ ^inner 



log orbital energy [J] 


Fig. 5.— Effect of relative orbital energies of the in¬ 
ner planets and the giant on the ejection probability. A: 

Chance of ejecting the giant rises as the orbital energy of the the 
giant planet relative to the inner planets falls. B: Where our exem¬ 
plar systems sit in the population of Kepler triple -transit systems 
( mass es estimated using a mass-radius relation, I Weiss & Marcvl 
120141) . For Kepler-18 and Kepler-58, we mark two values, corre¬ 
sponding to the masses estimated from the mass-radius relation 
(tail of arrow) and those corresponding to TTV/RV measurements 
(head of arrow), which were the ones used in the integrations. 

velocity (RV) surveys. A query of the Exop lanet Or¬ 
bit Database (http://exoolanets.org/. [Han et al.1 12014 
accessed 2015-05-16) revealed 354 K.V-detected planets 
with masses above 0.3Mj. Of these, 33 lie within 0.1 AU 
and 59 between 0.1 and 0.76 AU, 13 of the latter hav¬ 
ing e > 0.4. When weighted by their geometric transit 
probability, this sample has a higher fraction of hot to 
warm Jupiters than the KOI sample (see Table [I]), more 
in line with the ratio from our simulations. However, if 
we divide the warm Jupiters into two eccentricity bins 
at e = 0.4 (above whic h in-situ scattering is inefficient 
at exciting eccentricity (IPetrovich et al.l[2014f) . and be- 
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Fig. 6.— Orbital elements at the end of our integrations, for the Kepler-18 system with initial giant planet’s semimajor 
axis of 10 AU. A: Orbital elements of the giant planet when all inner planets are destroyed. Lighter symbols mark larger initial pericentre 
q. The red line marks the location of the giant planets at the beginning of the integrations. B: Orbital elements of surviving inner planets 
in systems where the giant was lost, marking multiplicity and merger state of the surviving systems. C: Evolving the distribution of giant 
planets from pan el A u nder tidal forces for 1 Gyr. Horizontal red lines mark the mean eccentricity attained through planet-planet scattering 
({Petrovich et all 12014 ). together with associated ler error bars assuming a Rayleigh distribution. We can produce highly-eccentric warm 
Jupiters at a < 1 AU with eccentricities higher than achievable through in-situ scattering of giant planets starting at these separations on 
circular orbits. D: Evolution of the surviving low-mass planets under tidal forces for 1 Gyr. The high eccentricities excited by interaction 
with the giant are unlikely to persist for the star’s main-sequence lifetime. E, F: As C and D but for 10 Gyr of tidal evolution. 
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low which tidal circularisation and/or interaction with 
the inner planets cannot reach), we find over 8 times as 
many hot Jupiters as eccentric warm Jupiters. This may 
point to a contribution from disc migration to the low- 
eccentricity warm Jupiter and hot Jupiter populations, 
although a detailed treatment of the differences between 
the RV and KOI samples is beyond the scope of this pa¬ 
per. 

When the incoming giant is ejected, the inner planets 
often experience some perturbation. Collisions of inner 
planets with each other or with the star are common, and 
the interaction with the giant often leaves systems with 
only one or two of the original three planets. The eccen¬ 
tricities of inner planets can be strongly excited (Fig [6)), 
and single survivors in particular can reach very high 
eccentricities. However, these eccentricities may not sur¬ 
vive in the long term, as tidal circularisation acts on 
the planets’ orbits on long timescales: Fig [6] shows that 
most eccentricities will decay to zero within 10 Gyr. In 
contrast, mutual inclinations of the inner planets are 
not strongly affected, although very flat systems do not 
retain their coplanarity: initially flat and moderately- 
inclined (few degrees) systems show similar inclination 
distributions after ejection of the giant (Fig [7]). 

Finally, we discuss the systems where the giant and 
at least one low-mass inner planet coexist at the end of 
the integration, shown in grey in Figure [3] In the over¬ 
whelming majority of these systems, the giant planet’s 
pericentre lies beyond the orbit of the outermost inner 
planet, and the final coexisting system looks much like 
the initial setup, retaining a highly-eccentric giant on 
a wide orbit with one or more low-mass planets close to 
the star; the number of inner planets may however be de¬ 
pleted by collisions. Interestingly, we do find a very few 
cases (23 in number) where at the end of the integration 
the giant planet’s semi-major axis is smaller than that 
of the outermost surviving inner planet. In 17 of these 
cases, all from the Kepler-58 simulations, both the giant 
and planet d he beyond 1 AU; in 6 of these, an additional 
planet remained at ~ 0.09 AU. In the remaining 6 cases, 
three each in the Kepler-58 simulations and the Kepler- 
18 simulations with the giant starting at 1.25 AU, the 
giant planet has collided with a b-c merger product, and 
the specific energy of the resulting body is sufficiently 
low that its orbit lies interior to that of planet d. Five 
of these systems survived integration for 10 Myr, and in 
all cases the giant planet’s pericentre is sufficiently small 
as to permit tidal circularisation and the formation of a 
hot Jupiter. In these five systems, the mutual inclina¬ 
tion is very high, oscillating around 90° and hampering 
the prospects for detection of both planets by transit 
photometry. However, these hot Jupiters with surviv¬ 
ing companions form only 0.5% of hot Jupiters formed 
in the Kepler-18, a = 1.25 AU integrations and 0.9% of 
those formed in the Kepler-58 integrations. Hence, while 
survival of companion planets is possible given the right 
conditions (viz. sufficiently massive inner planets), it oc¬ 
curs in only a tiny fraction of even these systems. We 
show the semi-major axes and eccentricities of planets in 
these coexisting systems in Figure@ highlighting the few 
systems in which the giant planet lies interior to one of 
the smaller ones. 


5. DISCUSSION 

5.1. The Evolutionary Context of our Simulations 

Our method assumes that the evolution of the system 
can be broken into three stages: an initial stage of excita¬ 
tion of the giant planet’s eccentricity; interaction of the 
highly-eccentric giant with any inner planets; and sub¬ 
sequent tidal circularisation of surviving planets’ orbits. 
We do not explicitly treat the initial phase, since the 
parameter space of mechanisms and perturbers is very 
large, and the combination of the small integrator step 
size needed to resolve the inner planets’ orbits conflicts 
with the long timesc ales (which can be over 10 s years, 
IWu & Lithwickll2011h needed to excite the eccentricity 
of the giant planet. The set-up for our simulations is 
probably most accurate for a scattering scenario, where 
the giant planet’s pericentre is impulsively changed to a 
very low value. In a Kozai or other secular scenario with 
a smoother eccentricity excitation, it is likely that the 
secular evolution will either continue until the timescale 
for secular evolution is comparable to the timescale for 
interaction between the inner planets and the giant (see 
FigHl), at which point our integrations begin; or that the 
interaction with the inner planets briefly halts the secu¬ 
lar cycle until they are destro yed (s imila r to the effects of 
general relativistic precession. IWu fc Murravf2003H . after 
which the secular cycle may resume. Note that destruc¬ 
tion of the inner planets can occur at pericentres wider 
than those at which the giant’s orbit actually overlaps 
the inner planets’ (Fig[3|. 

We also neglect any further effect of the body perturb¬ 
ing the giant planet during our integrations and after 
they have finished. This is again most accurate for a 
scattering scenario, where the swift reduction of the gi¬ 
ant’s apocentre following interaction with the inner plan¬ 
ets would decouple the giant planet from its original per- 
turber. In a secular or Kozai scenario, the secular cy¬ 
cles may resume after the inner system has been cleared, 
which will affect the statistics of hot and warm Jupiters 
we have estimated (Tables [3] and SJ) . Following the en¬ 
tire evolution of these systems from initial eccentricity 
forcing through to final tidal circularisation would be a 
fruitful avenue of future research. 

Although we have not treated the full evolution of these 
systems in this work, we can attempt to relate the out¬ 
comes of the integrations to the eccentricity excitation 
mechanism. In particular, a large semi-major axis of 
the giant planet increases the probability that it will be 
ejected instead of destroying the inner system. Driving a 
planet’s pericentre to very small distances by scattering 
from very wide orbits (note that to achieve a semi-major 
axis of 10 AU, the scattering event would h ave to tak e 
place at around 20 AU) is difficult dMustill et al. 1120III) , 
and the giant planets that we find vulnerable to ejection 
when they interact with the inner planets may be more 
likely to have been excited by Kozai perturbations from 
a wide binary companion. 

5.2. Robustness of our Findings 

The main result of our study—that giant planets with 
sufficient orbital eccentricity to become hot Jupiters de¬ 
stroy low-mass inner planets in the system—is robust to 
the masses of these inner planets, so long as they are 
not so massive as to eject the giant. In the absence of 
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Inclined Coplanar 



Fig. 7.— Distributions of mutual inclinations of surviving inner planets in Kepler-18 systems that lost their giant. A: 

With a spread of mutual inclinations up to 10°. We also mark the initial distribution (not to scale). B: Initially coplanar inner planets. 
In both cases the inner planets are not usually significantly excited, although there is a tail of inclinations above 10°. The distributions 
arising from the initially coplanar case and the initially inclined case are not significantly different. 


damping mechanisms that can separate and circularise 
orbits, the intersecting orbits of the giant and the inner 
planets lead to either collisions or ejections until orbits 
no longer intersect. In contrast, in very young systems, 
eccentricity can readily be damped by the protoplane¬ 
tary gas disc or by massive populations of planetesimals, 
helping to explain why Type II migration of giant planets 
does not totally suppress the formation of other planets 
in the inner parts of these systems: bodies thrown out by 
the gia nt ca n r ecir cularise and accrete outside its orbit 
(IMandell et alJ[2007f) . 

In our systems, in contrast to systems during the pro¬ 
toplanetary disc phase, gas is no longer present, and 
massive planetesimal populations are impossible to sus- 
tain close to the star for long time-scales dWvatt et all 
I2007D . Two additional sources of damping may play a 
role in these systems. First, debris may be generated 
in hypervelocity collisions between inner rocky planets, 
but integrations with the mass of the Kepler-339 plan¬ 
ets distributed among 100 smaller bodies did not show 
significant damping of the giant’s eccentricity. Second, 
tidal circularisation acts, but on timescales much longer 
than the time for planet-planet interactions to end in 
our systems. 

6. CONCLUSIONS 


We have shown that high-eccentricity migration of 
a giant planet to form a hot Jupiter necessarily leads 
to the removal of any pre-existing planets on orbits of 
a few tenths of an AU in the system, thus account¬ 
ing for the observed lack of close companions to hot 
Jupiters. This supports a high-eccentricity migration 
scenario for hot Jupiters, as migration through a proto¬ 
planetary g,as disc usually does no t fully suppress planet 
formation (Fogg fc Nelsonll2007allb]; Mandell et al.ll2007l : 


tormation (Fogg Nelson Zinu&.b: Man dell et al.ll/UUYt 

iKetchum et al.ll2011fc iQgihara et al l 20141) . We find that 


under high-eccentricity migration, when the giant’s peri- 
centre is sufficiently small to permit tidal circularisation, 
either the giant or the inner planets must be lost from 
the system. A very small fraction (< 1% even with 
favourable parameters) of the hot Jupiters we form do 
end up interior to a surviving low-mass planet, but this 
outcome is very uncommon: if such a low-mass close 
companion to a hot Jupiter were in future to be found, 
it would mean that in that system at least the migration 
almost certainly proceeded through a disc. When the 
giant planet does destroy the inner system, the interac¬ 
tion sometimes raises the pericentre of the eccentric giant 
planets sufficiently to prevent tidal circularisation, pro¬ 
viding a novel way of producing eccentric warm Jupiters; 
other giants whose pericentres are initially too high for 
tidal circularisation may be brought to populate the same 
region as they lose energy due to interaction with the in¬ 
ner planets. 

It is unknown which mechanism of eccentricity exci¬ 
tation dominates, be it scattering, the Kozai effect, or 
low-inclination secular interactions, but we expect that 
the inability of inner planets to survive in systems form¬ 
ing hot Jupiters will remain a robust result when future 
simulations coupling the evolution of the outer system, 
driving the giant’s eccentricity excitation, and the inner 
system are performed. 
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a [AU] 


a [AU] 


a [AU] 


Fig. 8.— Final (neglecting tidal effects) orbital elements of the giant planets (black) and low-mass inner planets (red) in the integrations 
in which both the giant and at least one of the inner planets are present after 1 Myr. “Inner” planets that now lie exterior to the giant are 
marked with large red stars, and planets in these systems are linked with black lines. 
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